Bose Condensates with 1/r Interatomic Attraction: Electromagnetically Induced 

"Gravity" 

D. 0'Dell (1) , S. Giovanazzi^, G. Kurizki«, V. M. AkuW 2 ) 
(1) Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel 
(2) Laboratoire Aime Cotton, CNRS II, Orsay Cedex 91405, France 
(February 1, 2008) 



OS 

On 
On 

o 

Q 



> : 
. 

On ■ 

o ■ 

<n : 

On ■ 
On ■ 

> ■ 

cd 



We show that particular configurations of intense off- 
resonant laser beams can give rise to an attractive 1/r inter- 
atomic potential between atoms located well within the laser 
wavelength. Such a "gravitational-like" interaction is shown 
to give stable Bose condensates that are self-bound (without 
an additional trap) with unique scaling properties and mea- 
surably distinct signatures. 

PACS: 03.75.Fi,34.20.Cf,34.80.Qb,04.40.-b 



In the atomic Bose-Einstein condensates (BECs) cre- 
ated thus far the atoms interact only at very short 
distance in good correspondence with the hard sphere 
model. The majority of the properties of these dilute 
gases can be understood by taking into account only 
two-body collisions which are characterized by the s-wave 
scattering length Q . A number of groups H have inves- 
tigated the fascinating possibility of changing the magni- 
tude and sign of the s-wave scattering length using exter- 
nal fields. The resulting condensates retain the essential 
hard-sphere, s-wave, nature of the interatomic interac- 
tion. 

Here we wish to introduce a qualitatively new regime 
of cold atoms in which the atom-atom interactions are 
attractive and have a very long range, varying as r" 1 [Q. 

We shall demonstrate that a stable BEC with at- 
tractive r~~ x interactions is achievable by irradiating the 
atoms with intense, off-resonant, electromagnetic fields. 
The atoms are then coupled via the dipoles that are in- 
duced by these external fields (in contrast to those in- 
duced by the random vacuum field responsible for the 
van der Waals-London interaction, which varies as r~ 6 , 
and leads to the usual hard sphere description) ||. 

Such an attractive potential can simulate gravity 
between quantum particles. Remarkably, this potential 
gives an interatomic attraction (depending on the laser 
intensity and wavelength) which can be some 17 orders 
of magnitude greater than their gravitational interaction 
at the same distance. 

This suggests it might be possible to study gravita- 
tional effects, normally only important on the stellar 
scale, in the laboratory. Particularly interesting is the 
possibility of experimentally emulating Boson stars ||: 
gravitationally bound condensed Boson configurations of 
finite volume, in which the zero-point kinetic energy bal- 
ances the gravitational attraction and thus stabilizes the 



system against collapse. 

In this letter we shall discuss the interplay of the usual 
hard-core interatomic potential with an electromagneti- 
cally induced "gravitational" one on a BEC using a vari- 
ational mean- field approximation (MFA). Two new phys- 
ical regimes with unique scaling properties emerge where 
the BEC is self-bound (no trap required). 

How can one realize the r _1 potential between neu- 
tral atoms? Consider the dipole-dipole interaction en- 
ergy induced by the presence of external electromagnetic 
radiation of intensity /, wavevector q, and polarization 
e. This energy can be written (in S.I. units) in terms of 
cartesian components i,j || 
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Here r is the interatomic axis, a(q) the isotropic, dy- 
namic, polarizability of the atoms at frequency cq, and 
Vij is the retarded dipole-dipole interaction tensor 



(5ij — 3f jf j) ( cos qr + qr sin qr) 
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where fj = Tijr. For a fixed orientation of the inter- 
atomic axis with respect to the external field, ([!]) and 
(^) give the well known r~ 3 variation of the interaction 
energy at near-zone separations (qr <C 1). The near zone 
limit of U (r) is strongly anisotropic. It was noted by 
Thirunamachandran jsj that when an average over all 
orientations of the interatomic axis with respect to the 
incident radiation direction is taken, the static dipolar 
part of the coupling (i.e., the instantaneous, non-retarded 
part r~ 3 (5ij — 3fjfjJ) vanishes. The remaining 'trans- 
verse' part is, in the near-zone, an attractive r -1 poten- 
tial. It is weaker by a factor of (qr) 2 than the r~ 3 term. 

However, thus far no scheme has been suggested 
wherein an average over all orientations is guaranteed 
for cold gases. We shall consider a spatial configuration 
of external fields which enforces the 'averaging out' of 
the r~ 3 interactions. A simple combination which en- 
sures the suppression of the r -3 interaction while re- 
taining the weaker r -1 attraction in the near-zone, uses 
three orthogonal circularly polarized laser beams point- 
ing along x,y,z ('a triad' — see Fig. 1). Let us mo- 
mentarily ignore interference between the three beams, 
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and only consider the sum of their intensities. In the 
near zone one can Taylor-expand Eqs. (|l|) and (||) in 
powers of the small quantity qr. Using the identity 
e* (±) (q)e^ ±) (q) = \ ((dy - q.q^) ± ie ijk q k ), with +(-) 
corresponding to left(right) circular polarizations, to- 
gether with Eqs. ([[]) and (||), the triad can be shown 
to give rise to the (near-zone) r _1 pair potential 
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Note that this interaction is attractive for any orientation 
(6, <fi) of r relative to the beams as long as the polariz- 
ability a{q) is real. 





FIG. 1. (a): Schematic depiction of a triad of lasers inci- 
dent upon an ensemble of atoms. This triad generates the 
attractive r" 1 potential given by Eq. (3), whose magnitude 
has the angular dependence shown in (b). 

If one wishes, the angular anisotropy in (^ can be can- 
celled to give a purely radial r _1 potential by combining 
a number of such triads with different orientations. It 
is convenient to define the orientation of each triad by 
the Euler angles (a,/3, 7) (?]], namely, a rotation of a 
about the z axis, followed by a rotation of f3 about the 
new y axis and finally a rotation of 7 about the final z 
axis. One configuration that cancels the anisotropy com- 
pletely uses 6 triads (18 laser beams) rotated through 
the following Euler angles: (0, 7r/4, 7r/8), (0, 7r/4, -7t/8), 
(0,7r/4,37r/8), (0,tt/4, -3tt/8), (0,0,tt/8), (0,0,-vr/8). 
The last two triads should be of half the intensity / of 
the others. Then the interatomic potential becomes 
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The main difficulty in realizing the near-zone r 1 po- 



tential is that the 



interaction survives due to the 



interference between different pairs of beams, whose 
contribution is proportional to the product of their re- 
spective field amplitudes. This difficulty can be over- 
come if one introduces frequency shifts between the laser 
beams. Spreading the frequencies w n of the E n laser 
fields (n — 1,2,3 for one triad or n = 1, 2 ... 18 for six 
triads) in intervals about the central frequency makes the 
r~ 3 interference terms in the interaction energy oc E n E^, , 



If these difference frequencies are much higher than the 
other relevant frequencies (e.g., collective oscillation fre- 
quencies) then the contribution of the interference terms 
to the mean field potential averages to zero. Typically 
these conditions hold for \u n — u) n >\ > 10 4 Hz. Angular 
misalignment errors, 5, between the orthogonal beams 
should satisfy 8 ^ qL (where L is the mean radius of 
the condensate) and intensity fluctuations should satisfy 
AI/I <C qL, in order to ensure the r~ 3 cancellation for 
the non-interfering terms oc J2 n \^n\ 2 - Although these 
oscillating r~ 3 terms do not contribute to the mean field 
potential, they can eject atoms from the condensate, but 
this process can be strongly reduced, as will be discussed 
at the end of this letter. 

A lower bound on the magnitude of the r _1 attraction 
is obtained by using the static rather than dynamic po- 
larizability, which for sodium atoms, say, has the value 
24.08 x 10~ 24 cm 3 . Thus, for a strongly off-resonant 
light intensity of I = 10 4 Watts/cm 2 one finds — u/r 
w —2 x 10 _19 eV, at r = lOOnm, the mean separation in 
a typical BEC. This is only around 10 of the magni- 
tude of the van der Waals-London dispersion energy at 
this distance. However in a system of many atoms the 
r _1 potential acts over the entire sample whereas the van 
der Waals-London interaction is only effective for nearest- 
neighbors and so the r _1 contribution to the total energy 
can become important. 

Our treatment of the many particle problem is based 
on a two-body potential V(r) — 4Ttah 2 5(r)/m — u/r, 
where the first term is the the pseudo-potential arising 
from the s-wave scattering (a is the s-wave scattering 
length and m the atomic mass). In order to write V 
in this form we require that the —u/r potential be suf- 
ficiently weak (compared with the mean kinetic energy 
per particle) so as not to affect the short-range hard- 
sphere scattering. This requirement certainly holds if 
a» \db a,, where Xdb is the de Broglie wave- 
length and a» = h?/mu is the Bohr radius associated 
with the gravitational-like coupling u. With the val- 
ues given above a* ~ 10 3 m, whilst for a typical BEC 
Ads ~ 10~ 5 — 10 -3 m and a ~ 3 nm. 

Consider now the application of this two-body poten- 
tial to a trapped dilute BEC gas well below the criti- 
cal temperature. We assume that the condensate ini- 
tially occupies a fraction of the wavelength of the laser 
so that the near-zone condition is valid (lasers in the far 
infrared, or microwave sources would satisfy this con- 
dition). The "zero-temperature" many-body problem 
leads, within the MFA, to the following equation for the 
order parameter \&(R, t) 
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(ri 7^ n')i oscillate at the difference frequencies \u) n — u) n >\. where V cx t(-R) = muj^R 2 /2 is for simplicity an isotropic 
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trap potential (which can be set to zero in certain cases — 
see below), and Vh(R) is the self-consistent Hartree po- 
tential 
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The order parameter ^(R, t) is normalized so that 
/ d 3 R |*(R, t)\ 2 = N, with N the total number of atoms. 
The usual Gross-Pitaevskii (GP) equation Q is recovered 
in the limit when u = 0. The MFA is valid when the sys- 
tem is dilute, i.e. pa 3 <C 1, with p the density. An 
additional condition on the MFA validity is that the r _1 
potential is weak, and this can be expressed as pa 3 ^> 1. 
This constraint, as in the related problem of the charged 
Boson gas [pj, means that many atoms must be present 
within the interaction volume a 3 . However, it is in fact 
the diluteness condition that turns out to be the stricter 
of the two, as one can check from the MFA results for 
the density (see Table 1 below). 

An analytical estimate for the mean radius of the N- 
atom condensate can be given using the following varia- 
tional wave function 



if A (i?) = N?(icX 2 lly* exp(-i?72A 2 ^) 
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where = yWpmujQ. The variational parameter A is 
proportional to the root mean square radius through 
\J (R 2 ) = \/3/2 XIq. This parameter is obtained by min- 
imizing the variational mean field energy 

H(X) 



N hojQ 
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where we have chosen the dimensionless u (proportional 
to the "gravity" strength u) and s (proportional to s-wave 
scattering length a) to be 
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The numerical factors are chosen to make the equation 
for A simple 



- A~ 4 + 1 + uX- 3 



sA- b = 



(10) 



This equation is equivalent to requiring that the vari- 
ational solution satisfies the following virial relation: 
-T + V cxt ~\U u ~lU s = 0, where T, V ext , U u and 
U s are the kinetic energy, the harmonic trap potential 
energy, and the internal energies due to the —u/r and 
hard-sphere interatomic potentials, respectively. This re- 
lation can be obtained from scaling considerations (see 
Ref. H for the case u = 0). 

The general asymptotic properties of the ground state 
solutions of Eqs. (||) and (||), as a function of (w, s), are 



summarized in the "phase diagram" of Fig. 2a for positive 
scattering lengths. In this diagram there are four asymp- 
totic regions: The non-interacting ideal region (/) and 
the ordinary Thomas-Fermi region (TF — O) are domi- 
nated by the balance of the external trap potential with, 
respectively, the kinetic energy and the repulsive s-wave 
scattering, and so are not sensitive to the —u/r poten- 
tial. The regions G and TF — G, which represent two 
new physical regimes for atomic BECs, are controlled by 
the balance of the gravity-like potential with either the 
kinetic energy (G) or the s-wave scattering {TF — G). 
Neither region is sensitive to the external trap, so that 
we can adiabatically turn it off (V ex t = 0) and access 
either the G or the TF — G region. 
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FIG. 2. (a): Contour plot of log(A), where A is the con- 
densate radius, in the parameter space log(u) versus log(s): 
darker shade corresponds to smaller A. The border separat- 
ing the TF — O and TF — G regions is given by s = u 5 ^ 3 and 
that separating the TF — G and G regions by s = u~ . (b): 
Mean energies per particle for large u (no external trap) as 
a function of the condensate radius. Curves are plotted for 
positive as well as negative values of the scattering strength 
s/u. For s/ii < —1/4 there is no minimum for a finite radius. 
The energy and radii units are u 2 Tiu>o and lo/u, respectively. 

Experimentally, direct signatures of the r _1 interaction 
come from the radius A and the release energy E re i — 
T + U s . The release energy is the kinetic energy that 
can be measured after the expansion occurring due to 
switching off the external trap and the laser fields [Q. 
Table 1 summarizes these quantities as well as the peak 
density p max in the four regions. 
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TABLE I. A comparison of the four asymptotic regions. 
The radius, A, and the release energy, E le \, are discussed in 
the text. p m ax is the peak density (at the center) of the con- 
densate. 
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We now focus on the properties of the two new regimes: 
a) In the TF — G region an analytic solution for the 
ground state of Eqs. (||) and (^J) is given by 



w 



2ir x 10 2 s - \ ql ~ 1, 



5, we obtain: 
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where Rq = ^Ja a*/2. Contrary to the ordinary Thomas- 
Fermi limit of the GP equation, the size of the condensate 
is fixed by the ratio of the coupling constants, AiraU 2 /m 
and u, and is independent of N. b) The G region, where 
only the r~ l attraction and kinetic energy play a role, is 
of particular interest since our system is then equivalent 
to a Boson-star (a system of gravitating Bosons) || in 
the non-relativistic regime. The mean field equations in 
this region are also identical to those describing a single 
particle moving in the gravitational field generated by its 
own wavefunction p0| . In both cases smooth bound so- 
lutions have been shown to exist |(],[l(| . This establishes 
the possibility of a stable self-bound (no external trap) 
r _1 condensate. 

The gravitational-like attraction does not induce "col- 
lapse" of the condensate, since, at short radii, it is al- 
ways weaker than the kinetic energy. This can be seen 
from the scaling of the kinetic energy (A -2 ) versus that 
of the r' 1 potential (A -1 ) in Eq. (§) and Fig. 2b. By 
contrast, this kind of instability can occur for negative 
scattering lengths when N exceeds a critical num- 
ber (N cr w 0.6 x |a|/Zo) because the mean energy due to 
scattering (A~ 3 ) is dominant at small radii. The u/r at- 
traction does reduce, when combined with the attractive 
scattering, the critical number to N cr sa 0.17 x y/a*/\a\ 
(see Fig. 2b for the critical case with su = —1/4). 

Finally, we estimate the losses of G or TF — G con- 
densates due to the r~ 3 oscillating interfering terms 
discussed above. Consider one of the possible oscillat- 
ing interfering terms A(r) cos(ili), where A(x,y,z) = 
—3 it ^2^5 and ^ is the difference in frequency between 
the two interfering lasers. Using Fermi's golden rule one 
can derive an expression for the rate of depletion of the 
condensate density I'J'I 2 due to creation of a pair of quasi- 
particles of opposite momenta (with k « mQ,/K) in 
the ideal homogeneous Bose gas 
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where |A(fc)| 2 is the angular average of the square of 
the Fourier transform of A(r) (0.1418u 2 /g 4 ). For our 
purposes it is sufficient to apply Eq. ( |l2|) at each point 
R (1*1 2 = |*| 2 (R)). Wc then find the following ap- 
proximations: dN /dt w u 5 s/TkJo/ (qlo) 4 in the G region 
and dNo/dt « u 7 / 2 S' 3 / 2 ^E^/(ql ) 4 in the TF-G re- 
gion. These expressions can be used to find conditions 
such that these loss rates are smaller than, say, the trap 
oscillation frequency ujq. Taking, e.g., £1 « 2n x 10 4 



for the G region d No/dt ~ 6 x 10 u>o, i.e. we need 
more than 10 5 atoms; for the TF region (e.g., s « 1), 
d N /dt « 6 x 10 3 w i i- e - we need more than 10 4 atoms. 



To conclude, the laser-induced attractive 



interac- 



tion can give rise to stable condensates with unique static 
properties. Their stability, long lifetime (low loss rates 
incurred by the r~ 3 oscillating terms) and lack of sensi- 
tivity to alignment errors or amplitude noise of the laser 
beams makes the experimental realization of such con- 
densates rather likely. Their fascinating analogy with 
gravitating quantum systems whose gravitational inter- 
action can be enormously enhanced by the field merits 
further research. 
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